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Abstract 

We propose a new interpretation of the recent result of the AMS-02 experiment, that Dark 



Matter (DM) can be two component because it shows flash damp at around 125 GeV of the 
positron energy. One of interpretations of that is primary source from annihilation into charged 
leptons of multi-component DMs. It is important to fix the flavor of final states in this case. 

This is achieved by imposing non-Abelian discrete symmetry as same as deriving the proper 

H ■ 

neutrino mixing. By assuming two gauge-singlet fermionic decaying DMs, we show that a model 

with non-Abelian discrete flavor symmetry, e.g. T\s, can successfully explain the damp of the 

positron fraction as well as the whole-range data. Few dimension six operators of universal 

leptonic decay of DMs are allowed in our model, since its decay operators are constrained by 

the T13 symmetry. We also show that the lepton masses and mixings are consistent with current 

experimental data, due to the flavor symmetry. 
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1 Introduction 

The latest experiment of Planck p] tells us that about 26.8 % of energy density of the universe 
consists of Dark Matter (DM). On the other hand, recent result of the indirect detection experiment 
of AMS-02 [2] is in favor of the previous experiments such as PAMELA [3] and Fermi-LAT [I] , which 
reported excess of positron fraction in the cosmic ray. Moreover, it smoothly extends the anomaly 
line of positron fraction with energies up to about 350 GeVLj. Along this line of thought, several 
papers have been released [6TPT3]. However, since AMS-02 as well as Fermi-LAT has also shown a 
flash damp of the line at around 125 GeV, that might imply an existence of lighter DM at this range 
if its damp has truly physical meaning 2 ]. 

In this paper, we show that this damp as well as the whole-range data can be explained by two- 
component fermionic decaying DMs. We introduce two pieces of decaying DMs (with mass of O(10 2 ) 
GeV and 0(1) TeV) into the framework of T 13 flavor symmetric scenario [15J. In our model, the 
flavor symmetry T 13 works (at least) in two ways: (i)It constrains interactions between DMs and the 
Standard Model(SM) particles. Since DMs (gauge singlet fermion X and X') couple with leptons 
by dimension six operators LELX^ ' /A 2 due to the T13 symmetry, these are leptophilic. DMs decay 
into leptons via these operators with the suppression factor A ~ 10 16 GeV, giving desired lifetime 
of DMs, T- 1 ~ ((TeV)VA 4 )- 1 ~ 10 26 sec HUE]. (ii)Flavor of the final states of DM decay is 
determined by the T13 symmetry. In this paper, we give a concrete example of the universal final 
states DM — > v e e + e~ / v [ Ji n + ji~ / v t t + t~ . 

The lighter DM is applied to explain the flash damp at around 125 GeV. The heavier one is, on 
the other hand, done to depict the main anomaly line up to 350 GeV that suggests the mass being 
larger than 700 GeV. Due to a specific selection rule by the flavor symmetry mentioned above, we 
show that two-component DM model is preferable for explanation of the AMS-02 result. In addition 
to that, we find a set of parameters that is consistent with the observed lepton masses and their 
mixings (especially rather large angle of #13) recently reported by several experiments [T8T - f2"^] . 

This paper is organized as follows. In the next section, we briefly mention the group theory of 
T13 and show the multiplication rules. After that, we construct mass matrices of the lepton sector 
in definite choice of T 13 assignment of the fields, and show that there exists a consistent set of 
parameters. In section 3, we show that desirable dimension six DM decay operators are allowed by 
T13 symmetry and that leptonic decay of the two DMs by those operators shows good agreement 



1 These observations can be, in general, explained by scattering and/or decay of the GeV/TeV-scale DM particles. 
Still leptophilic DM is preferable since PAMELA measured negative results for anti-proton excess [5] . 
2 This damp also occurs in the experiment of cosmic-ray proton and Helium spectra [14] . 
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Table 1: Characters of Ti 3 . ^ is defined as the complex conjugate of £j 



with the cosmic-ray anomaly experiments. Section 4 is devoted to the conclusions. 



2 Lepton masses and mixings with T13 flavor symmetry 



First of all, we briefly review the non-Abelian discrete group T13, which is isomorphic to Z13 x Z% 
[T5 | 1251427] . The T13 group is a subgroup of SU(3), and known as the minimal non-Abelian discrete 
group having two complex triplets as the irreducible representations. We denote the generators of 
Zi 3 and Z 3 by a and b, respectively. They satisfy a 13 = 1 and ab = ba 9 . Using them, all of T 13 
elements are written as g = b m a n , with m — 0, 1, 2 and n — 0, 1, • • • ,12. 
The generators, a and 6, are represented e.g. as 

/o 1 o\ / 



V 



1 

1 



/ 



p 



\ 



V 



p 3 








where p = e 



2i7r/13 



(2.1) 



Q 
P ) 



These elements are classified into seven conjugacy classes, and T13 has three singlets 1& with 
k = 0, 1, 2 and two complex triplets 3i and 32 as irreducible representations. The characters are 
shown in Table (TJ where £1 = p + p 3 + p 9 , £ 2 = p 2 + p 5 + p 6 , and u = e 2l7r//3 . The multiplication rules 



of the T13 group are shown as 



3i®3i = 


= 3i©3i©3 2 , 


3i©3i 


= 3iffi3i©3 2 , 


3 2 <8> 3 2 = 


= 3 2 ©3 2 ©3i, 


3 2 ©3 2 


= 3 2 ©3 2 ©3i, 


3i®3 2 = 


= 3 2 ©3 2 ©3i, 


3i©3 2 


= 3i©3 2 ©3i, 


3 2 ® 3i = 


= 3i©3i©3 2 , 


3i©3 2 


= 3i©3 2 ©3 2 , 


3i®3i = 


= ^ l fc ©3 2 

fc=0,l,2 


©3 2 , 




3 2 ©3 2 = 


= ^ l fc ©3i 


©3i. 





(2.2) 

fc=0,l,2 

The tensor products between singlets are found as lo®lo = li<8>l2 — lo, li©li = I2, I2C3I2 = li- 
The tensor products between triplets and singlets do not have an affect on their property. 

Let us move on discussion the lepton masses and mixings in the setup shown in Table HJ Here, 
Q, U, D, L, E, H(H') and X(X r ) denote left-handed quarks, right-handed up-type quarks, right- 
handed down-type quarks, left-handed leptons, right-handed charged leptons, Higgs bosons, and 
gauge singlet fermions, respectivelyj. Here notice that X and X' is Majorana- and Dirac-type DM, 
respectively. Due to the T 13 flavor symmetry in addition to an appropriate choice of the additional 
Z 3 symmetry, triplet Higgs bosons -H"(3i) and i?(3 2 ) couple only to leptons, while T13 singlet Higgs 
bosons -ff'(lo,i,2) couple only to quarks. Therefore, the mass matrices of quark sector are not 
constrained, while those of lepton sector are determined by the T 13 symmetry. For the neutrino 
sector, since the Yukawa couplings LHX^ are forbidden by the T 13 symmetry, the left-handed 
Majorana neutrino mass terms are derived from dimension five operators LHLH. Here notice that 
X(X') has dimension six operators LELX^ ) and mass terms. For the matter content and the T13 
assignment given in Table HI the charged-lepton and neutrino masses are generated from the T 13 
invariant operators 



'Y 



V2a e ELH c (3 2 ) + \F2b e ELH\Z x ) 
+ ^-LH(3 2 )LH(3 2 ) + ^LF(3 2 )LF(3!) 
+ jLH^LH^ + h.c, (2.3) 

where H c = eH* and the fundamental scale A = 10 11 GeV (For the lifetime of DM, A/V\ ~ 10 16 
GeV is required to obtain the desired lifetime, where A is coupling constant of DM decay operators). 



3 All the assignment and particle contents are the same as our previous work [15] except the DM sector. 
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Table 2: The T13 and Z3 charge assignment of the SM fields and the Majorana DM X and the Dirac DM 
X' where u = e 2i7r / 3 . 



After the electroweak symmetry breaking, the Lagrangian Eq. fl2.3j) gives rise to mass matrices of 
charged leptons M e and neutrinos M v 
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where the vacuum expectation values (VEVs) of the Higgs bosons are defined as 
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Now we give a numerical example. By the following choice of parameters 



Vl 

a e 

b„ 



4.269 GeV, v 2 = 161.1 GeV, v 3 = 78.62 GeV, 
10 GeV, v 2 = 168.2 GeV, v 3 = 0.04836 GeV, 
0.01057, b e = 0, a u = -8.220 x 10~ 5 , 
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the mass matrices Eqs. ( 12. 4ft and ( 12.51) give rise to mass eigenvalues and related observables as 
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0.511 MeV, m„ = 105.7 MeV, 
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1777 MeV, m ul = 6.324 x 10" J eV, 
1.078 x 10~ 2 eV, m u3 = 5.046 x 10~ 2 eV, 
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2.43 x 10~ 3 eV 2 , 



2.83 x 10" 4 eV, ^m vi = 5.49 x 10~ 2 eV, 



(2.8) 



and the mixing matrices 
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which are all consistent with the present experimental data |28l|29]. In particular in the case of 
U e L = 1, the mass matrices Eqs. (12 .4p and (12.51) require normal hierarchy m u i < m U 2 < vn v3 of 
the neutrino masses and £/Jf 7^ 0. A comprehensive analyses of the T 13 symmetric models have 
been done by several authors [30]. In the following analysis, we assume these parameters to discuss 
decaying DMs. 



3 Decaying dark matter in the T13 model 

It is well known that the cosmic-ray anomalies measured by AMS-02 [2], PAMELA [3] and Fermi- 
LAT [4J can be explained by DM decay with lifetime of T" 1 ~ 10 26 sec. If the DM (X and X' in 
our case) decays into leptons by dimension six operators LELX^ ) /A 2 and A ~ 10 16 GeV, such long 
lifetime can be achieved. In general, however, there exist several gauge invariant decay operators 
of lower dimensions: dimension four operators inducing too rapid DM decay, and dimension six 
operators including quarks, Higgs and gauge bosons in the final states, which must be forbidden 
in a successful model. By the field assignment of Table [2j most of all the decay operators listed 



Dimensions DM decay operators 

4 LH C X^ 

5 
6 LELXV, WHLH c X { -'\ {H c ) t D ll H c EYX { !\ 

QDLX { '\ UQLX('\ LDQX('\ Uj lt DE'fX^, 

D^H C D^LX('\ D^D^H C LX('\ 

B llu La^H c X^'\ W^La^r a H c X^ 

Table 3: The higher dimensional operators which cause decay of X and X' up to dimension six |31j . 
Bfiu, W® u , and D^ are the field strength tensors of hypercharge gauge boson, weak gauge boson, and the 
electroweak covariant derivative. 

in Table |3] [31] except for LELX^ ' are forbidden due to the T13 symmetry Q With the notation 
U = (vi,£i) = {U eL ) ia (u a ,£ a ) and E { = {U eR )ipEp (i = 1,2,3, a, ft = e,fi,r), the four-Fermi decay 
interaction is explicitly written as 

A; ^ 



Alccay = -^ J2^ E ^ X +( X ^ X ') + h - C - 
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= ^ E ( w2(i_1) ) E ^)*<* (^ (^ 



x [(P a P fl ^) (^PrX) - (l 7 P^) (p a P R X)] +(X^ X') + h.c, (3.1) 

where in the case of X' decay, the factor (u; 2 ^ 1 -*) has to be multiplied because of the multiplication 
rule of the T 13 flavor symmetry. As seen from Eq. (13. ip . decay mode of the DM particles X(X') 
depends on the mixing matrices U e (i,R), which are given in Eq. ( 12. 9p . 

Next, we consider the branching fraction of the DM decay through the T13 invariant Lagrangian 
Eq. ( 13. ip . Due to the particular generation structure, the DM particles X(X') decay into several 
tri-leptons final state with the mixing-dependent rate. The decay width of DM X per each flavor is 
defined as 

r> 7 = t(x -> uj+q) + t(x -)■ j^ep-), (3.2) 



4 Notice that WHLH C X^ > and H^HXX cannot be forbidden. Here we assume these couplings of these surviving 
terms to be enough tiny. However, we expect that these terms can be forbidden by modifying the other flavor 
symmetries with higher order elements such as £(81) of T,(3N 3 ) [32], A5 [33], the double covering of Ac, [17], and Tig 
ofTjv. 



and the decay width r a/ 3 7 is calculated as 
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The decay width of X' named T' a o is obtained by replacing X — y X '. Note that the factor (w 2 ^ -1 - 1 ) 
is needed in the case of X' decay. The differential decay width is required to get energy distribution 
function dN/dE e ±, and is expressed as 

d^a/3-y _ I Ax I Tn x 



-x 2 ((6 - 2x)U aM + (15 - Ux) U ai p\ (3.5) 



dx 48(4vr) 3 A 4 ' 

where x = 2E e +/m x . Here we have omitted the masses of charged leptons in the final states. In 
both X and X' DM cases, the flavor dependent factor V a ^ gives a factor three if one takes the sum 
of flavor indices a, f3 and 7. Therefore, the branching fraction of each decay mode is given by 

Br(X -»- vJ+£-) + Br(X -)■ Vjfyq) = ~ (C/^ 7 + C/ a7/3 ) . (3.6) 

The DM mass m^ and the total decay width T x = J2 a b 7 ^a/3-y ar e chosen to be free parameters in 
the following analysis. In the case of X', we can obtain it by replacing X — y X' . 

Given the differential decay width and the branching ratios, the primary source term of the 
positron (electron) from DM decay at the position r of the halo associated with our galaxy is expressed 
as 

Q(E,r) = n x (r)T x J2MX^f)(^f) + (X -y X') , (3.7) 

where [dN e ±/dE]f is the energy distribution of positrons (electrons) from the decay of single DM with 
the final state /, and E is the energy of e ± . We use the PYTHIA [34] to evaluate the distribution func- 
tion [dN e ±/dE]f. The non-relativistic DM number density n x (r) is written by n x (r) = px{ r )/ m x 
with the DM profile p(r). In this work, we adopt the Navarro- Frank- White (NFW) profile [35J, 

t \ r &( r e + r c) 2 ,„ Q x 

pNFw{r) = p Q — 1—. vj-, (3.8) 

r(r + r c ) A 

where p Q ~ 0.40 GeV/cm 3 is the local halo density around the solar system, r is the distance from 
the galactic center whose special values r ~ 8.5 kpc and r c ~ 20 kpc are the distance to the solar 
system and the core radius of the profile, respectively. 
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Figure 1: The positron fraction [2 J and [3] predicted in the leptonically-decaying DM scenario with X13 
symmetry. We fix the lighter DM mass to be mx =244 GeV and the heavier DM mass to be mx' =1360 
GeV. As for the DM decay width used in the fit, see the text. 

As seen from Eqs. (I2.9P and (J3.3J) . the DM decays into r^ as well as e ± and /x ± in the equal 
rate. As a result, pure leptonic decays give dominant contributions, and it is consistent with no 
anti-proton excess of the PAMELA results [5] . We follow Ref . [16] for diffusion model describing the 
propagation of positrons and electrons [3BT438] , and backgrounds [361I39] . 

Result for AMS-02 The positron fraction is depicted in Figure Ufor the scenario of the leptonically 
decaying DM with T13 symmetry. We can obtain the best fit point in terms of four parameters mx, 
m x> ', Tx, Tx' by x 2_n t and the values are 



m x 
m X ' 



244 GeV, 
1360 GeV, 



r -i 

1 x 



V2S 



r 



X' 



1.2 x 10 z " s, 
2.4 x 10 27 s, 



(3.9) 
(3.10) 



where we used 35 data of AMS-02 from higher energy. We obtain the value of x 2 — 14.18 (31 d.o.f) 
at the best fit point. In this figure, we fix to the acquired masses and decay widths, assuming that 
the number density between X and X' is the same. The result is shown with the experimental data 
of AMS-02 and PAMELA. One can see a small damp around E rj 125 GeV. This is the result 
of mult i- component DMs and the fixed flavor of final states by T13 symmetry. A sharper damp is 



8 



expected if we have larger branching ratio for directly produced positron. 

4 Conclusions 

We revisited a decaying DM model with a non-Abelian discrete symmetry T 13 , and extended it to 
the two-component DM scenario by adding an extra DM X' in order to explain the flash damp of the 
anomaly line at around 125 GeV. We have shown that all the observed masses and mixings in the 
lepton sector are derived from the T 13 model consistently Due to also the specific selection rule of 
T 13 , we have found that both of DMs have the universal decay coming from dimension six operators 
that gives a promising model in current indirect detection searches. In this paper, we have explicitly 
shown a numerical examples of a consistent set of parameters in the DM mass, their lifetime, and 
mass matrices of the lepton sector. We have found that the experimental small damp of AMS-02 
can be explained by two-component DM scenario of nix =244 GeV with T^ 1 = 1.2 x 10 28 s and 
nix' =1360 GeV with T x ) = 2.4 x 10 27 s, assuming that X and X' have equal number density. 

Acknowledgments 

We would like to thank Prof. Shigeki Matsumoto for a crucial advice. T.T. acknowledges support 
from the European ITN project (FP7-PEOPLE-2011-ITN, PITN-GA-2011-289442-INVISIBLES). 
The numerical calculations were carried out on SR16000 at YITP in Kyoto University. 

References 



[2] 
[3] 

[4] 

[5] 

[6] 

[7] 



P. A. R. Ade et al. [ Planck Collaboration] , |ar2£ivil303 i 5076| [astro-ph.CO]. 

M. Aguilar et al, Phys. Rev. Lett. 110, 141102 (2013). 

O. Adriani et al, Nature 458 (2009) 607. 

M. Ackermann et al. [Fermi LAT Collaboration], Phys. Rev. Lett. 108, 011103 (2012) 



jarXiv: 1109. 052 1 [astro-ph.HE]]. 

O. Adriani et a/., Phys. Rev. Lett. 102 (2009) 051101. 

P. -F. Yin, Z. -H. Yu, Q. Yuan and X. -J. Bi, larXiv:1304.4"T28l [astro-ph.HE]. 



Q. Yuan and X. -J. Bi, arXiv: 1304.2687 [astro-ph.HE]. 
J. Kopp, larXiv: 1304. 1T841 [hep-ph]. 

9 



[9] A. De Simone, A. Riotto and W. Xue. larXiv: 1304. 13361 [hep-ph]. 

[10] Q. Yuan, X. -J. Bi, G. -M. Chen, Y. -Q. Guo, S. -J. Lin and X. Zhang, I arXiv: 1304 .14821 [astro- 
ph.HE]. 

[11] M. Ibe, S. Iwamoto, S. Matsumoto, T. Moroi and N. Yokozaki. [arXiv: 1304. 14831 [hep-ph]. 

[12] T. Linden and S. Profumo. larXiv:1304.179~Il [astro-ph.HE] . 

[13] H. -B. Jin, Y. -L. Wu and Y. -F. Zhou, larXiv: 1304. 19971 [hep-ph]. 

[14] O. Adriani et al. [PAMELA Collaboration], Science 332, 69 (2011) [ arXiv: 1103.40551 [astro- 
ph.HE]]. 

[15] Y. Kajiyama and H. Okada, Nucl. Phys. B 848, 303 (2011) [arXiv:1011.5753l [hep-ph]]. 

[16] N. Haba, Y. Kajiyama, S. Matsumoto, H. Okada and K. Yoshioka, Phys. Lett. B 695, 476 
(2011) |arXiv: 1008.47771 [hep-ph]]. 

[17] K. Hashimoto and H. Okada. larXiv:1110.3640l [hep-ph]. 

[18] T2K Collaboration, K. Abe et al, Phys.Rev.Lett. 107, 041801 (2011), [1106.2822]. 

[19] J. K. Ahn et al. [RENO Collaboration], Phys. Rev. Lett. 108, 191802 (2012) jarXiv: 1204.06261 

[hep-ex]]. 

[20] DAYA-BAY Collaboration, F. An et al, Phys.Rev.Lett. 108, 171803 (2012), [1203.1669]. 

[21] DOUBLE-CHOOZ Collaboration, Y. Abe et al, Phys.Rev.Lett. 108, 131801 (2012), 
[1112.6353]. 



[22] CHOOZ Collaboration, M. Apollonio et al, Eur. Phys. J. C27, 331 (2003), p^e x/0301017]. 

[23] Palo Verde Collaboration, F. Boehm et al, Phys. Rev. D64, 112001 (2001), |hep-ex/01070~0~9~] . 

[24] MINOS Collaboration, P. Adamson et al, Phys.Rev.Lett. 107, 181802 (2011), [1108.0015]. 

[25] W. M. Fairbairn and T. Fulton, J. Math. Phys. 23 (1982) 1747. 

[26] S. F. King and C. Luhn, JHEP 0910(2009) 093. 

[27] For a review of non-Abelian discrete symmetry, H. Ishimori, T. Kobayashi, H. Ohki, H. Okada, 
Y. Shimizu and M. Tanimoto, Lect. Notes Phys. 858, 1 (2012); H. Ishimori, T. Kobayashi, 
H. Ohki, H. Okada, Y. Shimizu and M. Tanimoto, Prog. Theor. Phys. Suppl. 183 (2010) 1. 

[28] J. Beringer et al [Particle Data Group Collaboration], Phys. Rev. D 86, 010001 (2012). 

[29] D. V. Forero, M. Tortola and J. W. F. Valle, Phys. Rev. D 86, 073012 (2012) [ arXiv:1205.4"0T8l 
[hep-ph]]. 

10 



[30] K. M. Parattu and A. Wingerter, Phys. Rev. D 84, 013011 (2011) [arXiv:1012.2842l [hep-ph]]: 
G. -J. Ding, Nucl. Phys. B 853, 635 (2011) [arXiv:1105.5879 [hep-ph]]; C. Hartmann and A. Zee, 
Nucl. Phys. B 853, 105 (2011) [ arXiv:1106.0"333l [hep-ph]]; C. Hartmann, Phys. Rev. D 85, 
013012 (2012) |arXiv:1109.5T43l [hep-ph]]. 

[31] F. del Aguila, S. Bar-Shalom, A. Soni and J. Wudka, Phys. Lett. B 670, 399 (2009) 



[32 
[33 

[34; 

[35 
[36 
[37 
[38 
[39 



[arXiv:0806.0876 [hep-ph]]. 



H. Ishimori and T. Kobayashi, Phys. Rev. D 85, 125004 (2012) [arXiv: 1201 .34291 [hep-ph]]. 

F. Feruglio and A. Paris, JHEP 1103, 101 (2011) |arXiv:1101.0393l [hep-ph]]. 

T. Sjostrand, S. Mrenna and P.Z. Skands, JHEP 0605 (2006) 026; Comput. Phys. Commun. 
178 (2008) 852. 

J.F. Navarro, C.S. Frenk and S.D.M. White, Astrophys. J. 490 (1997) 493. 

E.A. Baltz and J. Edsjo, Phys. Rev. D 59 (1998) 023511. 

D. Hooper and J. Silk, Phys. Rev. D 71 (2005) 083503. 

D. Maurin, F. Donato, R. Taillet and P. Salati, Astrophys. J. 555 (2001) 585. 

C. Pallis, Nucl. Phys. B 831 (2010) 217. 



11 



